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DEFINITION OF SYMBOLS 


EMR 

emr 

FM 

f 

g 

h 

LM 

i 

MB 

R 


u c/o 

W(r) 

Wi(i) 

w 2 (f) 

M 


*1 

2 


°EMR 

*R 


Random variable, stage effective mixture ratio 
Outcome of EMR 
Nominal mainstage fuel 
fcR F 

g:EMR -*■ Rj^ 
h:EMR - R F 
Nominal mainstage LOX 
i € Rp 
Fuel bias 

Random variable, outage (R^U Rp) 

Set of all possible fuel outage 
Set of all possible LOX outage 
Outcome of R 

Nominal stage mixture ratio at cutoff 
Nominal tanking ratio 

FM 

Probability density function of R 

Frequency of occurrence of I 

Frequency of occurrence of f 

Mean value of payload potential, excluding 
the effects of outage 

Mean outage 

Normal probability density function 

Standard deviation of payload potential, excluding 
the effects of outage 

Standard deviation of EMR 

Standard deviation of outage 
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SECTION I. INTRODUCTION 


Outage is defined as the amount of propellant (fuel or oxidizer, 
but never both) which cannot be burned because of an insufficient 
quantity of the other propellant. Outage obviously reduces the vehicles' 
capability potential in the form of unwanted inert weight at cutoff. An 
outage will occur when a variation exists between the targeted and actual 
tanked mixture ratio and/or between the nominal and the actual burned 
mixture ratio. 

Since outage is clearly a random phenomenon, any payload study 
involving outage will be statistical in nature. The tanked mixture ratio 
and burned mixture ratio are assumed normal with known mean and 
standard deviation. The effective mixture ratio (EMR) is the statistical 
combination (root sum square) of these two ratios. In Section II the 
probability density function (pdf) of outage is developed analytically as 
a function of LM, FM, “eMR^c/o ant * bias. The equations show 
how the fuel bias helps shape the pdf of outage and hence, how it influ- 
ences the mean and standard deviation of outage. 
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SECTION II. THE PROBABILITY DENSITY FUNCTION OF OUTAGE 


Consider the maps h:EMR-*Rp and g:EMR-*RL defined as follows 
(Figure 1): 



FIGURE 1. STAGE MIXTURE RATIO 


and 


f = h(emr) 


= FM + MB 


LM 

emr 


emr > 


LM 

FM + MB 


l = g(emr) 

= (^H/o| fm+mb 


LM 

emr 


emr 


< LM 
~ FM + MB 
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By c on structi on > the functions g and h are 1-1 and they map the 
disjoint sets 

LM . LM 

emr < ... . . . emr > . . 

FM + MB FM + MB 

onto the disjoint sets 

' ^F* 


respectively. 

Hence. 

W,(r) 


the irequency of LOX outage occurrence may be given by 
= |5£-! • r « R^. zero elsewhere 


where q is the probability density function (pdf) of EMR and 



r 


is the absolute value of the Jaccbian of the inverse of the transformation g 
(References 1 and 2). 


Similarly for fuel 


W 2 (r) = q [h _ 1 (r)J | , r e Rp . zero elsewhere 

Now 

r c R = R ^ U Rp 

and 

R l 0 R F = <j> 

implies that 

P(r c R) = P(r c R L ) + P(r e R f ) 

which is to say that the pdf of R is 
W(r) = W,(r)+ W 2 (r) 


= Tifg-Mr)] 


dgl 1 

dr 


t rifh’ 1 (r)] 


dh 


-l 


dr 


r c R. zero elsewhere. 
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Let r « R^> then there exists an 


emr s 


FM + MB 


such that 


r = g(emr) 


“c/o |(FM + MB ) - - 5 -LK: 


= - LM (FM + MB )emr . 

ILp U«p 


Solving for g“* (r) 


g' 1 (r) = emr 


u c/o LM - u T r 


7 - r — ~ . r c R t , zero elsewhere 

u c/0 (FM + MB ) L 


Jn” ^ ~ Up* 

-* — = » rt R, . zero elsewhere. 

dr r u c y c (FM + MB) *-> 

Similarly 

h " 1(r> = FM + ~ iAB - ~ r ’ r * R F * zero el8ewhere 


dh- 1 _ [ 

I dr _ " f 


LM 

(FM + MB - r) 2 


» r c Rp . zero elsewhere. 


Hence. 


W(r) = iltg" 1 ^)] 


u (FM + MB) 


+ T1 [h - 1 (r)] 


LM 

(FM + MB - r) 2 


( u , LM - u T r \ 

u 1(FM 7 b H 


tr EMR^ r 


+ exp I- 

<r EMR^r L 


z^emr” 


LM 1 

.FM + MB - r ~ U T; 


Z^EMR* 
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which is the form of the pdf of outage [W{r)l used for digital computation 
of the mean outage n 

or = ) rW(r> 


and standard deviation of outage 


[J r 2 W(r)J- u R 2 


li 


in the computer prc gram that is discussed in the following sections. 
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SECTION in. COMPUTER PROGRAM LISTING 

*** I NP'jT ************* ************* 

*S I GE= 0 NE STANDARD DFVI AT ION OF EFFECTIVE. MIXTURE RATIO * 

* nrOaCTAf.P CUTOFF MlXT'jRF RATIO * 

*pm=MA I NST AGE FURL LOAD. * 

*L f< *M A I NST AOP LOX LOAD. * 

»viQ_pijpL p J AS * 

* * * OUTPUT *********************** * * 

*o = pp.ST OUL r S (LR*') * 

♦DP-ppo^AP I L I TY OF P or LFSS OCCURRING. 

*FP=FRFQU c NfY OF OCCURRANCF OF R. 

PCSIOUAL. * 

*Sir,;< = STAf\OARn OF V I A T I ON OF RFSIDUAL. * 

#do, jv- ruTOFF LIMIT For PR * 

* I PR T = OPTION FOR DETAIL PRINT OF R.FR.PR I F ( I PRT .F-F . 1 ) * 

♦NoaSoSs i Tmpu * 

******************************** 

p f a l L v » v R 

rov*iQA| /cy*y /LPI *NNY 
oimpmSION A ( 2 5 0 0 » 4 ) 


o I MFNSlON TOX(12.3).IOY(12*3 )»LappL( 7,3) »LHFA0( M ♦ j SYM{3 ) ,L( * ) 
o i *<*PNS I ON XL ( a l ,XR ( * ) * YP ( > ) ,YT ( * I *TOZ (1 2 ) »L° ( “> ) 

NAMCLIST / NaMpI/ STOF*UcO»FM*LM*Mq,UT»CI*r3*PRLTN’*lPRT 
N aMfL I ST / .N aI , 'f 2 / T OX . TOY »L ARPL .LHFAO, I OZ » I SYrf , XL , yR ♦ YR * YT , L , LP 
* ♦ tPLOT *Nr aSfs 

P F A O ( S » N A M F y ) 

WD tTF( 6*NA mc 2 ) 

T F ( I PLOT ,LT, 1) GO TO *> 

FAIL 0 AMR AV ( 3 SR ) 
v MYs? 

*** SYSTEM SUBROUTINE **** 
r A! L SFOllTV 

FALL PITF2V I80.5S0.900»9 A .l , 36»1 »LHFAO*NfRR) 

c rf>MTl M HP 

WP T TF ( ft * R I 0 ) 

0° 200 IJ = 1 * N C A S F S 
Rfao( 5»NAMFl ) 

UT*LM/FM 
PP = 0. 

UP,0. 

RRFRcO , 

R 3 P P a 0 t 
PPPFPsO, 
c IOR?sO. 

DRa S . 

■»*n, 

PPPrVaO. 

A ( 1 ♦’ )a 0 , 

A ( 1 t’lsO. 

A ( 1 * * ) = 0 , 
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H =1 , / (MGF*tSQRT( ?.M%141ft9?7) ) ) 
r'»aUT/(Uro*(FM+MB> ) 

VRTTc-(6,NAMPl ) 

r>0 too isi,?*>Cn 
I T = I 
’•* = T - 1 

PRsR - RPR?V 
PPPFVsP 

r? s | ( ( (UrO*LM)-(iJT*R) ) / ( 'JCO* { FM+Ml * ) ) l"L'T )**? / ( *> .*< I np*« ! fie ) 

r&N=(LM/(FM+\ie-R) )-UT 

rts r4N*r4N/(?.*SinF*siGF) 

re = A"S(LfV( ( pM+MR-R ) **? . ) ) 

rp = ri*(cxP(-r?l*r's + fxP ( -cu ) *nft 1 

T F ( ) no TO 10 

PP=po+OR*( (CR+A(y*?1 )/?.) 

UR = |JR + PR*(R*CFR+A <M*2) )/?. ) 

P?FR=R*R#FR 

RRFRsRRFR + DR*( (R?FR+PR?FR) /?. ) 

pr?fr=r?fr 
r I oRpsRRfR - (UR*UR) 

10 O0MT1NUF 
a ( T *i )*R 

A ( 1 »? )=PP 
A ( I * 0 ) =PR 

TP (PR ,r,T. PPL I V ) PO TO 10] 

TOO rONTlNUF 
1 01 00NTINUF 

ft TORsSOPT ( SK,R9) 

WR T TP ( ft » R 0 1 ) UR » ft I GR »P* 

]P( IPRT ,LT. 1 ) no TO 140 
PO 10? T = 1 * IT •* 

10? WR T TP ( 6 * *5 04 ) T . A ( I * 1 ) * A ( I ♦ ? ) » A ( ! « K ) 

140 ]P( IPLOT.I.T, 1) GO TO 1 50 

nO ?0 J*i ♦*» 

JJ*J+i 
LPi=LP( J1 

r #** SYSTPM SUPROUTTNF **** 

C ALL OUIKOVtL ( J) * ISYM( J) « I DX ( 1 » J ) » I OY ( 1 ♦ J ) • I I » A (1 * 1 ) , A ( 1 ♦ J j ) ) 
PALL R!T^2V(250»990,900»9‘ ,36tJ *LHRAO*N pR. > 

CALL R I TF2V( 210,960, 980,90,1,42,1 *LAPFL<1»J) ,NpRR> 
f CO^TINUF 
1*0 CONTINL'R 

VRTTP(ft»F?0) 

'•’R r T p ( ft » A ? 1 ) STGF 
VR]TP(ft,ft??)UrO 
wo T 1 p( ft , ft ?? ) pm 
, .'.'PTTP(ft,R?4)LM 
»'R lTP(ft,ft?ft)^ft 
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r?TTc-ns*e?M 

■>'RTTrt5* e ?7)Sir,R 

'/."5TTc((C,,r->b )UP 
S'PTTFCfetRO? ) 

>(,Q roNl I N JE 

I F ( ! PLOT .IT. 1) STOP 
'*'R T T F ( i 4 , * l'J ) 

TALL R ITr2V( 89,550 ,900,9? ,1 ,36 ♦l.LHFAn.NrRR) 
r »** 5YSTFM SUBROUTimF **** 

r ALL CLEAN 
crop 

501 FORMAT ( 4H UR= . E16 .8 . 1 5h SORT < SI GR2 ) * » E 16. 8 , 6H PR*.F16.8! 

502 PCPVAK ////12H FND OF EASE) 

5?3 FORMATt 105X,5H E4N=,F14.7) 

504 format (V; *I4*iSF14.7) 

51 0 FORMAT ( 1 Hi ) 

520 fORMaTUhI *////*?7X*5MTNr At/) 

571 FoovaTUBH STANDARD DEVIATION OF EFFECTIVE MIXTURE RATIO =,ro.4> 
522 FORMAT ( 26H STAGE CUTOFF MIXTURE RA T lO = ,2f)X,R9.4 ) 

573 F0R*'AT(?2H MAINSTAGF FUEL LOAO= *2 IX »F14.4 ) 

524 FOP y AT ( 22H ‘*A ! NSTAGF LOX LOAD =*21X»F14.4) 

525 FORMAT ( 12H FUEL BIAS = . 31X-.F14.4 ) 

525 FORMAT ( //// ,1H * 27X *6H0UTPUT ♦ / ) 

527 FORMAT ( 32H STANDARD DEVIATION OF R* SIDUAL= »9X t f 14. 7 ) 

47* FdpmaT<75H EXPECTED VALUE OF RESIDUAl ( MEAN > * *6X tFlA. - ’ 1 

FND 
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SECTION IV. COMPUTER PROGRAM SAMPLE RESULTS 


Table I shows the results of nine computer runs that were made 
to support a design phase. Space Shuttle sizing study. Total run time 
for the nine cases was 2. 3 min on the Univac 1108 computer. 


TABLE I. PROGRAM RESULTS 



Input Data 

Output Data 

Case 

°EMR 

U c/o 

LM 

FM 

MB 

**R 

or 

R 

1 

0.021 

mm 

1. 399. 175 

233, 196 

300 

1413 

1543 

psf! 

0. 021 

m 

1, 399, 175 

233. 196 

o 

o 

' J 

1224 

1264 


0.021 

B 

1, 399, 175 

233, 196 

800 

1243 

1046 


0.021 


1. 399, 175 

233, 196 

900 

1255 

974 

5 

0.021 


1, 399, 175 

233, 196 

1000 

1281 

915 

6 

0.021 


1, 399. 175 

233, 196 

1100 

1320 

869 

B 

0.021 

B 

1. 399, 175 

233, 196 

1200 

1377 

836 

IS 

0.021 

m 

1. 399, 175 

233. 196 

1400 

1498 

799 

o 

0. 021 
1 

6. 0 

1. 399, 175 

233, 196 

1600 

1654 

7P8 


The computer program has the option of generating plots of the 
probability density function and distribution function of outage. Computer 
plots corresponding to case two (Figures 2 and 3) and to case six 
(Figures 4 and 5) are shown. 
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Frequency of Occurrence, W(r) 



FIGURE 2 . PROBABILITY DENSITY FUNCTION, W(r)--CASE TWO 
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W(r) 



FIGURE 3. DISTRIBUTION FUNCTION— CASE TWO 
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Frequency of Occurrence, W(r) 



FIGURE 4. PROBABILITY DENSITY FUNCTION, W(r)--CASE SIX 
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W(r) 



Outage (lbm) 


FIGURE 5. DISTRIBUTION FUNCTION- -CASE SIX 
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SECTION V. CHOOSING THE OPTIMUM FUEL BIAS 


Let M and Z be the mean and the standard deviation of payload 
potential, excluding the effects of outage. M and Z are the statistical 
combinations (usually root sum square) of payload potential variations 
caused by the randomness of thrust, ISP, propellant loads, and aero- 
dynamic forces. 

By an application of the Central Limit Theorem (Reference 3), 
the total standard deviation, including outage effects on payload 
potential, is closely approximated by 



Within a given sigma probability, k; the total vehicle payload 
potential is 

M - [|ij^ ±k/(r R 2 + Z 2 ] 

Now, given that LM, FM, ^EMR* u c/o are known, ^ and ar^ 
are uniquely determined for a given fuel bias (Section II). Thus 
optimizing the fuel bias amounts to choosing the fuel bias that 
minimizes the payload potential loss 

Hr + k/r R » ♦ S* 

within a given sigma probability, k. 

Table II shows the payload potential loss associated with each 
of the output data of the sample computer run (Table I). The 3Z used 
was 4151 lbm. The optimum fuel bias for this study was determined to 
be 1100 lbm (Table II). 
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TABLE II. PAYLOAD POTENTIAL LOSS (lbm) 
VERSUS FUEL BIAS 


Case 

MB 

+ / (3 V* + 4151 * 

1 

300 

7630 

2 

600 

6846 

3 

800 

6441 

4 

900 

6332 

5 

1000 

6258 

6 

1100 

6222 


1200 

6227 


1400 

6291 


1600 

6431 


It is worth noting that the optimum fuel bias is not highly 
sensitive to reasonable values of 2. Several fuel bias optimization 
studies were made for this Space Shuttle configuration with 32 ranging 
from 1000 to 6000 lbm. The resulting optimum fuel biases were between 
1000 and 1300 lbm. However, the optimum fuel bias is sensitive to the 
total propellants tanked (LM + FM) and to EMR. 
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SECTION VI. SUMMARY 


This paper presents an accurate and efficient method of deter- 
mining the optimum fuel bias for a bipropellant liquid rocket. Basically, 
the paper shows that the mean and standard deviation of outage are 
uniquely determined for a given fuel bias (Section II), and how, because 
of this, probable loss in payload resulting from outage can be minimized 
by the proper choice of a fuel bias. 
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